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Abstract 

Static flexoelectric effect in a finite sample of a solid is addressed in terms of phenomenological 
theory for the case of a thin plate subjected to bending. It has been shown that despite an explicit 
asymmetry inherent to the bulk constitutive electromechanical equations which take into account 
the flexoelectric coupling, the electromechanical response for a finite sample is " symmetric" . " Sym- 
metric" means that if a sensor and an actuator are made of a fiexoelectric element, performance 
of such devices can be characterized by the same effective piezoelectric coefficient. This behavior 
is consistent with the thermodynamic arguments offered earlier, being in conflict with the cur- 
rent point of view on the matter in literature. This result was obtained using standard mechanical 
boundary conditions valid for the case where the polarization vanishes at the surface. It was shown 
that, for the case where there is the polarization is nonzero at the surface, the aforementioned sym- 
metry of electromechanical response may be violated if standard mechanical boundary conditions 
are used, leading to a conflict with the thermodynamic arguments. It was argued that this conflict 
may be resolved when using modified mechanical boundary conditions. It was also shown that the 
contribution of surface piezoelectricity to the fiexoelectric response of a finite sample is expected 
to be comparable to that of the static bulk contribution (including the material with high values of 
the dielectric constant) and to scale as the bulk value of the dielectric constant (similar to the bulk 
contribution). This finding implies that if the experimentally measured fiexoelectric coefficient 
scales as the dielectric constant of the material, this does not imply that the measured fiexoelectric 
response is controlled by the static bulk contribution to the fiexoelectric effect. 

PACS numbers: 77.22.-d, 77.65.-j, 77.90.+k 
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I. INTRODUCTION 



The flexoelectric effect consists of a linear response of the dielectric polarization to a 
strain gradient. This is a high-order electromechanical which is expected, in general, to be 
rather weak. However, some of its features make this effect to be of interest from both funda- 
mental and applied points of view. This has stimulated recent intensive experimental^^-and 
theoretical-^^activity in the field. 

On the fundamental side, it is of interest that this effect cannot be considered just as a 
non-local generalization of the piezoelectric effect. In contrast to the later, the flexoelectric 
effect (response) is controlled by 4 mechanisms of different physical nature, the contributions 
of which can be comparable^. 

On the practical side, of first importance is that this effect, in contrast to the piezoelectric 
effect, is allowed in centro-symmetric material. It is believed that it is the flexoelectric effect 
that is responsible for the generation of an electric field in acoustic shock waves propagating 
in centro-symmetric solids^. It was recently shown that this effect plays also an essential 
role in electromechanical properties of materials with a moderate level of electronic and 
ionic conductivity^. However, the most applied interest is focused on the "piezoelectric 
metamaterial" - composites made of non-piezoelectric components, which exhibit effective 
piezoelectric response generated due to the flexoelectric effect. The work in this direction was 
initiated by pioneering experimental studies by Professor Cross with coworkers- and was 
later also supported by theory^. Presently, (Ba,Sr)Ti03-based composites have been shown 
to yield effective piezoelectric coefficients comparable to those of commercial piezoelectric 
ceramics^. It was argued, based on the constitutive equation for the flexoelectric response, 
that a mechanical sensor made of such metamaterials should exhibit a very unusual property. 
Specifically, in contrast to piezoelectric based devices, it will not behave as an actuator-i^^. 
There are several reasons to question such statement. First, already in the 60 's of the past 
century, the group of Professor Bursian reported experimental data on BaTiOs crystals^^ and 
gave arguments based on equilibrium thermodynamics^, which contradict this statement. 
Second, the existence of a linear sensor-not-actuator may come into conflict with the general 
principles of thermodynamics. 

The goal of this paper is to address theoretically this conflict situation to demonstrate 
that despite an explicit asymmetry of the constitutive equations for the bulk flexoelectric 



2 



effect, when this effect is characterized in a reahstic finite sample the apparent asymmetry 
of the electromechanical response will vanish. In particular, this implies that the aforemen- 
tioned piezoelectric metamaterial should exhibit the identical piezoelectric constants when 
characterized in "direct" and "converse" regimes. We will demonstrate this for two leading 
contributions to the static flexoelectric response: the contribution of static bulk ferroelec- 
tricity and that of surface piezoelectricity. 



II. STATIC BULK FLEXOELECTRICITY 



In the present paper, being interested in the static or quasi-static situation, we are fac- 
ing three contributions to the flexoelectric response, which are associated with (i) static 
bulk ferroelectricity, (ii) surface piezoelectricity, and (iii) surface fiexoelectricity''. In this 
section we will discuss the problem of the relation between the direct and converse effects 
for contribution (i), reserving the next section for contribution (ii). In view of vanishing 
practical importance of the surface flexoelectric effect (it is expected not to be enhanced on 
high-dielectric-constant materials which are the only ones suitable for applications) we will 
not be treating it in this paper. 

The static bulk flexoelectric effect is customarily described by the following free energy 
(density) expansion (see e.g.-): 

771 r> T-) fijkl ( Yi dUij dPk\ Cijkl 
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where Pi and Uij are the polarization vector and strain tensor, respectively, and where the 
Einstein summation convention is adopted. We will consider this thermodynamic potential 
as having the differential dF = EidPi + aijduij. Then calculating the electric field, Ei, and 
stress tensor, aij, as variational derivatives of the free energy of the sample, given by the 
integral of F over its volume, one arrives at the following constitutive equations: 

Ek = XkjPj - f^i^^^^ (2) 

and 

o'ij = fijki-Q^ + CijkiUki- (3) 

The first equation describes a linear polarization response to strain gradient (direct flexo- 
electric effect). The second one describes the converse flexoelectric effect, implying that to 



get a "mechanical yield", spatial inhomogeneity of the polarization is needed. From this, one 
might infer (as customarily done in relevant papers) that the application of a homogeneous 
electric field to a sample will not lead to its deformation. Even being nearly evident, in 
reality, the last statement is not correct. 
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FIG. 1: Plate of the material exposed to bending and the reference frame used in calculations. 

Let us show this for the flexural mode. Consider, a (001) plate of a cubic material 
of thickness h in the reference frame specified in FiglH with the X and Y dimensions 
being L and h respectively. To make the analysis transparent, we allow only a cylindrical 
bending of the plate about OX2 axis. To simplify the discussion further, we set, for the 
moment, C1122 = C1133 = 0. In such simplified model, the plate bending is associated 
with dun/dxs 7^ 0, whereas U22 = M33 = so that Eq.(j2]) suggests the appearance of P3 
component of the polarization controlled by /1133 component of the flexoelectric tensor. To 
address the reversibility of this effect, one should check if the application of an electric field 
normal to the plate will cause its bending. A straightforward way do this is to derive the 
equation of balance of the bending moment for the plate^'*' subjected to a homogeneous 
electric field E normal to its suface by integrating Eq.(l3]) across a YZ cross-section of the 
sample: 



At mechanical equilibrium, the Ihs term must be equal to the minus the component of the 
mechanical moment of the external forces applied to the lefthand (with respect to the cross- 
section of the integration) part of the plate,— M2. Without the first rhs term, this equation 
describes that bending of the sample caused by this moment. To identify the role of this 
term, we first evaluate it using integration by parts: 
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where (P3) is the averaged polarization induced by the field in the bulk of the plate. In 
doing so we assume that the polarization changes continuously from its bulk value to zero on 
the plate boundary. If, however, one explicitly considers nonzero polarization at the sample 
surface (as was done previously in Ref-), then one should revise the traditional boundary 
conditions of the elasticity theory—. We will return to this issue later in the paper. Since the 
spatial scale of the polarization variation at the interface is much smaller than the thickness 
of the plate, with a good accuracy (P3) ^ P, where P - polarization in the bulk. Thus, the 
equation for the moment balance can be rewritten as 

rh/2 

- M2/6 + finshP = ciiii / uuzdz. (6) 

J-h/2 

It is clear from this equation that the application of a homogeneous electric field to the 
plate is equivalent to that of an external bending moment. Thus, we conclude that a finite, 
mechanically free (M2 = 0) sample, placed in a homogeneous electric field, will be bent. 
This conclusion is closely related to that drawn by Eliseev et al^. These authors have shown 
that a ferroelectric plate with the out-of-plane orientation of the spontaneous polarization 
should exhibit spontaneous bending due to the fiexoelectric coupling. It was found that this 
effect is controlled by a factor j'^l^^^{dP3/dz)zdz which was calculated using a numerical 
solution for the polarization profiles P{z) in the sample. Here it is also worth mentioning 
that the bending effect addressed, though being proportional to a component of the bulk 
fiexoelectric coefficient and the bulk value of the induced polarization, is actually controlled 
by forces applied to the surface of the plate. 

It is instructive to illustrate quantitatively the "symmetry" of the direct and converse 
fiexoelectric effects in a finite sample for a situation which is readily mathematically track- 
able. We will consider the case of a (001) plate of a cubic material in symmetrical fiexural 
mode with the polarization P normal to the plate and homogenous in its bulk. In the case 
of symmetric bending, the curvature of the plate in all crossections normal to it, G, is the 
same. Enjoying the results of the theory of thin plates^ we can express the components of 
the strain tensor in terms of this curvature: 

un = U22 = zG] U33 = -z — G] U12 = U23 = ui3 = 0. (7) 

Cll 

Integrating the free energy density, Eq.(IT]), with strain coming from Eq.([7]), over the plate 
thickness, one finds the free energy density per unit area of the plate as a function of P and 
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G: 

$ = ^hP' + ^G' - 2hPGifn33 - —/nil). (8) 

Z Z Cll 

JJ ^ '^11 + CIICI2 - 2C^2 /g-j 

6 Cll 

where Dg is a coefficient controlling the flexural rigidity of the plate for this kind of bending. 
In derivation of Eq.([H]), we have again used Eq.(|5]) (this gave a factor of 2 in the coupling 
term from this equation). Similar expression for the free energy of flexoelectric plate in the 
cylindrical bending mode was offered by Bursian and Trunov^^, based on purely symmetry 
arguments. Minimizing Eq. ([8]) with respect to P and G one finds the equation for the direct 
and converse effects for the plate in symmetric flexural mode: 

G = — /ipiE. (10) 

D = P = 2/ipiG. (11) 
where the electric displacement, D, and 

cii/1133 — C12/1111 ,^ . 

/^pl = X33 (12) 

Cll 

can be treated as an effective flexoelectric coefficient of the plate. The flexural response, 
given by Eq. (|TOl) . is compatible with the results obtained by Eliseev et al^ for the case of 
spontaneous bending of thin plates with the blocking boundary condition for the polariza- 
tion. 

Obviously, elements of such plate will work as both actuators and sensors. If round 
pieces of the plate with central loading and symmetric free-edge side support are used as 
elements of a piezoelectric metamaterial, electromechanical properties of the latter will be 
characterized by a single effective piezoelectric coefficient ^33. Using the relation between 
the cross-section curvature, G, and the maximal deflection, ^max, for symmetric bending of 
a circular plate: 

GR^ , > 

Smax ^ • K^'^J 

(where R is the radius of the plate) and using Eq. f lTOj) one readily finds 

43 = (14) 
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III. CONTRIBUTION OF SURFACE PIEZOELECTRICITY 



As was recognized at the first thorough treatment of the flexoelectric response^, the 
polarization response to a strain gradient in a finite sample, generally speaking, may not be 
fully controlled by the contribution of the bulk static fiexoelectricity, even in materials with 
high values of the dielectric constant (high-i^ materials). The competing effect that is due 
to surface piezoelectricity, was not, however, properly addressed theoretically. Thus, it is 
not clear if it can in fact compete in high-K materials with the static bulk fiexoelectricity. 
In high-i^' materials, the contribution of the static bulk fiexoelectricity is enhanced, since 
it scales as the dielectric susceptibility (cf. Eg. (ITT]) and Eq.( lT^ ). At the same time, the 
effect associated with surface piezoelectricity originates from the presence of the interface 
adjacent layers where the piezoelectricity is induced by the inversion-symmetry-breaking 
effect of the interface. Since the sign of the effective piezoelectric coefficients of the layers 
on the opposite sides of a plate should be opposite (as controlled by the orientation of the 
surface normal), bending of the plate should result in dipole moments in these layers, the 
sign of which are the same. The dipole moment in a layer is proportional to the strain in it, 
which, in turn, is proportional to the product of the strain gradient and the plate thickness. 
Having calculated the resulting change of the average polarization of the whole system, this 
will give rise to a net polarization proportional to the strain gradient. From this reasoning 
it is not obvious that such response will be enhanced once that the dielectric constant of the 
bulk of the material is high. However, such reasoning does not provide a proper vision of 
the whole effect. In what follows, we will show that such enhancement does take place and 
the considered bulk and surface contributions to the flexoelectric response can be readily 
comparable in high-i^' materials. We will also address the problem of the relation between 
the direct and converse effect for the mechanism related to the surface piezoelectricity. 

To be specific we will address these problems for the case of symmetric bending of a thin 
(001) plate of a cubic material. We will model the effect by considering a system consisting 
of a plate of an "ideally homogenous" material (i.e. its material parameters are the same 
throughout the plate) with the bulk flexoelectric effect being neglected and two thin surface 
piezoelectric layers (Fig 12]). The thickness of each layer. A, is much smaller that that of the 
plate, h. The top layer is characterized by piezoelectric moduli /isss and /isn = h^22i whereas 
for the bottom layer these moduli have the same the absolute value but are of opposite sign. 
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FIG. 2: Model for the contribution of surface piezoelectricity to the flexoelectric response of a non- 
piezoelectric material. The surface layers of thickness A model the surface adjacent (atomically 
thin) layers of the material where the piezoelectricity is induced by the symmetry breaking impact 
of the surface. 

We also ascribe to these layers an out-of- plane component of the dielectric constant, ^a- 
Let us find the extra free energy associated with the top piezoelectric layer when the plate 

is symmetrically bent with a cross-sectional curvature, G, and when out-of-plain polarization 

in the layer equals P\. We assume no difference in the elastic properties of the piezoelectric 

layers and the plate. 

We start with the free energy density in the layer defined as 

Fa = ^Pl- Px [/i333M33 + ^311 + M22)] + ^ (^n + M22 + ^33) + C12 (M11M33 + M22M11 + M22M33) 

(15) 

where Cn = Cim, C12 = C1122, and a is the inverse dielectric susceptibility of the layer if it 
were fully mechanically clamped. Because the plate is thin we set in the layer 

Mil = ^22 = hG/2. (16) 

As for M33, we find it from the condition that the surface of the film is mechanically free, 
dFx/du33 = 0: 

C12 , ^ , ^333 n r-tr,-. 

U33 = hC H Pa- (17) 

Cll Cll 

Inserting Eqs.f ll6p and f|T7|) ) into Eq. ffTSj) and multiplying the result with A, one finds the 
free energy of top piezoelectric layer. The energy of the bottom piezoelectric layer is the 
same since this layer differs from the top layer by the sign of piezo-moduli and by that of 
the strain; in the expression for energy these signs cancel each other out. Thus, for the 
contribution of the two piezoelectric layers to the free energy of the system, we find: 



$x = 2A 



Z Cll 



$1 (18) 



where xa = (« — ^333/cii)~^is the true (under the mixed mechanical conditions) dielectric 
susceptibility of the piezoelectric layers and $i is their mechanical bending energy. 

To describe the direct flexoelectric response we use the equation of state for the polar- 
ization in the piezoelectric layer, d^x/dPx = Ex {Ex is the electric filed in the layer), the 
condition of continuity of the electric displacement, D, in the layer, and the short-circuit 
condition. This leads to the following set of equations: 

Px = XxEx + ehG (19) 

D = efEf = eoEx + Px (20) 

2\Ex + hEf = (21) 

where Ej, ej, Eq are the electric field in the bulk of the plate, its dielectric constant, and 
the dielectric constant of the free space, respectively, and 

e = Xx{h3u-—h,3,). (22) 

Cll 

Solving this set of equations we find the relation for the direct flexoelectric response: 

D = 2eXG. (23) 

where 



hef 



(24) 



2Xef + hex 
and ex = eo + x\- 

To describe the converse flexoelectric response we present the elastic energy of the system 
at fixed Px in the form 

$ = - 2hXeGEx. (25) 

When writing this equation we have neglected the elastic energy of the surface layer, $i, 
compared to that of the plate and have simplified flTI?]) down to Px ~ X\Ex- This ap- 
proximation means that we neglect feedback effect of G on Px- This effect will yield some 
renormalization of -D^, but practically such a renormalization is negligible indeed. Minimiz- 
ing $ with respect to G and applying the electrostatic relation used above we arrive at the 
set of equations 

^ 2/iAe , , 

G = -^E, (26) 
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(27) 



2\Ex + hEf = E{h + 2A) 



(28) 



where E is the apphed (average) electric field. This set leads us to the equation for the 
converse fiexoelectric effect in the system: 



When writing this relation, only the leading terms to within a small parameter X/h were 
kept. 

The following remarks are to be made concerning the results obtained. First, the relations 
obtained for the contribution of the surface piezoelectricity into the fiexoelectric response, 
Eqs. (!23|) and (!29|) . are identical to those obtained for the case of static bulk fiexoelectricity, 
Eqs. ( fTTl) and ( fTOl) . to within the replacement /ipi =^ e\. Thus, all conclusions about the 
direct-converse-effect symmetry drawn in the previous section for the static bulk fiexoelec- 
tricity hold for the contribution associated with the surface piezoelectricity. Second, in a 
high-K material the latter contribution scales as its dielectric constant (similar to the case 
of the static bulk fiexoelectricity). Formally, this follows from the expression for e, Eq. (!24|) . 
Taking into account that the thickness of the surface piezoelectric layer is expected to be of 
the order of the lattice constant, for realistic values of the plate thickness, e can be evaluated 
as 



Note that there is no reason to consider ex in a high-K material (typically it is a "regular" 
of incipient ferroelectric in the paraelectric phase) to be high, since the special interplay of 
the atomic forces responsible for the high value of the bulk permittivity will be inevitably 
destroyed in the surface layer. Such enhancement of a surface driven effect by a factor of 
the bulk permittivity looks surprising. However, the physical mechanisms behind this effect 
can be identified. 

For the converse effect it is quite transparent. The bending of the system is controlled 
by the value of the field in the surface layer (cf. Eq. fl26|) ). Due to its small thickness, this 
field is enhanced by a factor of ej/ex , compared to the applied field. 

For the direct effect, the explanation is less straightforward. This time, the bending 
creates polarization in the surface layer. Because of the inhomogeneity of the system, the 




(29) 



e = 




(30) 
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short-circuiting does not guarantee the absence of the electric field in it so that the polar- 
ization in the surface layer induces a depolarizing field both in itself and in the bulk of the 
plate. It occurs that if the surface layer is thin enough whereas £f/sx is large, the polar- 
ization response is controlled by the depolarizing field in the bulk of the plate. This way 
the polarization response of the system becomes sensitive to the bulk value of the dielectric 
constant. 

Another important conclusion is that, taking into account the aforementioned effect of 
enhancement, one expects both contributions to the flexoelectric response discussed to be 
of the same order of magnitude even in high-K materials. These contributions would be 
comparable if Xe/sf were about the typical value of the components of the fiexoelectric 
tensor /jj^/ , 1 — 10 V (see e.g.-). For the "atomic values" of the entering parameters 
(A = 0.4 nm and e = 1 C/m^) and sx/eq = 10, we evaluate Xe/sf ^ 4 V to find that this is, 
in fact, the case. 

IV. DISCUSSION AND CONCLUSIONS 

The analysis presented clearly demonstrates that starting from continuous constitutive 
electromechanical equations one can derive the relations for the direct and converse fiexo- 
electric effects in a finite sample, which exhibit the symmetry required by thermodynamics. 
On the practical level such symmetry implies that a piezoelectric meta-material based on 
the fiexoelectric effect will exhibit the same effective piezoelectric coefficient in the testing 
regimes for the direct and converse piezoelectric effect. This is in conflict to the belief of 
those dealing with such meta-materials^ii^. 

At the same time, it is also clear that the analysis, based on the constitutive electrome- 
chanical equations presented above, is valid only for the situation where the polarization at 
the plate surface can be treated as continuously changing from its bulk value to zero. If it 
is not the case, formally following this analysis one readily flnds that the aforementioned 
symmetry is violated. For example, for the case of free boundary conditions for polarization 
{dP/dz = at the boundary), corresponding to nonzero polarization at the sample surface, 
the equation of mechanical equilibrium, Eq. (jlj), implies the absence of the converse effect. 
Such conclusion would be fully consistent with that by Eliseev et al^ who argued that the 
manifestation of the converse flexoelectric effect in a plate is strongly dependent on the 
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boundary conditions for the polarization. At the same time, there is no reason to expect 
that the free boundary conditions for the polarization will suppress the direct fiexoelectric 
effect. Thus, if we followed the calculating scheme employed by Eliseev et al^ (and used 
in SectJTT]) we would find, for the free polarization boundary conditions, the absence of the 
converse fiexoelectric effect in the presence of the direct effect. This would make an apparent 
contradiction between the results obtained from the continuous constitutive equations and 
those obtained from thermodynamics. 

We suggest the following resolution to this contradiction. The point is that incorporating 
the fiexoelectric coupling into the free energy density of a material leads to a modification 
of the boundary conditions for the bulk constitutive equations. Eliseev et al^ have derived 
modified boundary conditions for the polarization, however these authors have postulated 
that the classical mechanical boundary conditions are not affected by such incorporation. 
However, as was recently shown by one of the authora^^, generally, the mechanical boundary 
conditions should be modified as well. It has been shown that such boundary conditions 
reduce down to the classical mechanical boundary conditions when the polarization vanishes 
at the surface. This justifies the calculations based on the classical mechanical boundary 
conditions, which we have presented in SectUll For the general case, the problem of the 
converse fiexoelectric effects should be revisited with the correct mechanical boundary con- 
ditions which contain the surface value of the polarization. We expect that such treatment 
will yield the results consistent with the symmetry between direct and converse fiexoelectric 
effects in a finite sample dictated by thermodynamic arguments. 

Another important conclusion follows from the results obtained in Sect JIIII There, it was 
shown that the contribution of the surface piezoelectricity to the fiexoelectric response is 
expected to be comparable to that of the static bulk contribution (including the material 
with high values of the dielectric constant) and to scale with the bulk value of the dielectric 
constant (similar to the bulk contribution). The latter statement actually implies identical 
(or at least similar) temperature dependences of these contributions. Note that in earlier 
publications^ it was hypothesized that these depedences are expected to be different. Based 
on this hypothesis, the fact that the experimentally measured fiexoelectric coefficient scales 
as the dielectric constant of the material was customarily taken as an indication that the 
measured fiexoelectric response is controlled by the static bulk contribution to the fiexoelec- 
tric effect. The results from Sect JIIII essentially change the situation. Now one can state 
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that the fact that the experimentally measured flexoelectric coefficient scales as the dielectric 
constant of the material does not imply that the measured flexoelectric response is controlled 
by the static bulk contribution to the flexoelectric effect. 
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